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Studies of non-interacting lattice fermions give an estimate of the size of discretization errors and finite size 
effects for more interesting problems like finite temperature QCD. We present a calculation of the thermodynamic 
equation of state for free domain wall and overlap fermions. 



1. Introduction 

The thermodynamics of fermionic systems are 
certainly affected by lattice regularization in a 
finite volume. In numerical calculations of the 
QCD equation of state at currently accessible 
lattice spacings, the data tend to overshoot the 
continuum ideal gas limit. This general pattern 
is consistent with long established results [Q for 
non-interacting lattice fermions. This suggests 
that for any new lattice fermion formulation, one 
should first study the thermodynamics of the non- 
interacting system to gain insight on the potential 
cutoff effects to be expected when the interactions 
are turned on. For a recent discussion and an 
example of this approach, see the review talk of 
Karsch from LATTICE '99 §. 

2. Energy density of free fermions 

The energy density of a relativistic gas of non- 
interacting free fermions is given by the integral 

esBim) = — / -p dp 1 

Jo 1 + exp{^yp^ + rn^/T) 

which is known analytically in the massless limit: 
£sb(0) = (77rV60)r4. This is often called the 
Stefan-Boltzmann (SB) limit from the depen- 
dence. One important feature of this calculation 
is that a continuous set of momenta are summed 
to compute the energy density. However, sim- 
ilar calculations for free fermions on the lattice 
are summed only over the finite set of discrete 
lattice momenta. So, before we judge the rela- 
tive merits of various free lattice actions based 
on comparisons with the SB limit, it is useful to 



first understand consequences of restricting free 
fermions with a continuum dispersion relation to 
the allowed lattice momenta. 




N4 = (aT)-^ 

Figure 1. Discretized energy density of free 
fermions with am = 0, 0.1, • ■ • , 1. 

As this is a free field calculation, the scale is 
set explicitly by specifying the dimensionful tem- 
perature T. The infrared cutoff given by the 
spatial extent of the lattice can be specified in 
terms of the dimensionless aspect ratio: ^ = VT^ . 
The ultraviolet cutoff given by the lattice spac- 
ing can by specified in terms of the dimension- 
less A'4 = (aT)"-'^. If we choose a,V so that the 
product CiV4 is an integer, then the discretized 
analogue of (||) is 



~ + cxp(aiV4 VpM^) 



10 o o o o ^ 

o 

0.8 - 

CD 

-.^ 0.6 - 

T 
T 

^ 0.4 - 

J o 

0.2 - 



2 4 6 8 

^ (aspect ratio) 

Figure 2. Discretized massless energy density vs. 
aspect ratio. 

where the sum is over —n/a < p < n/a with 
spacing Ap = 2TT/CN4a. 

In figure |l|, we display the both the SB (|l|) and 
discretized (Q) energy densities at fixed aspect ra- 
tio C = 4 for various A'4 and masses from am = 
to 1 spaced in 0.1 intervals. For high temper- 
atures (small A'4), the mass of the fermions be- 
comes irrelevant and the SB (solid) lines converge 
on £sb(0). In this same limit, the number of al- 
lowed lattice momenta becomes small so the dis- 
cretized (dashed) lines are cutoff and converge on 
zero. For this aspect ratio, the SB and discretized 
values differ by more than a percent for iV4 < 4. 

In figure |l|, the massless discretized energy 
density falls short of the SB limit, approaching 
0.989 as iV4 00. For massless fermions, even 
small momenta contribute at all temperatures, 
and these are cutoff by the aspect ratio C = 4. 
In figure |^, we show the massless discretized en- 
ergy density vs. aspect ratio in the A'4 00 limit. 
So, from just the continuum free dispersion rela- 
tion and a simple discretization procedure, lattice 
simulations should not expect to reproduce the 
continuum energy density to better than a few 
percent on any lattice smaller than 16"^ x 4. 

3. Comparing free Wilson fermions 

The energy density of free Wilson fermions has 
been known for a long time |l[] and its large over- 
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Figure 3. Energy density of free Wilson 

fermions. 

estimation of the SB limit at smaller values of A'4 
has often been used as an argument against the 
applicability of Wilson thermodynamics at cur- 
rently accessible lattice spacings. Because stan- 
dard domain wall and overlap fermions use Wil- 
son terms to give mass to the doubler species, we 
might anticipate that similar effects may be seen 
in the energy density of these fermion formula- 
tions. So, we would like to first understand free 
Wilson thermodynamics. 

We show the energy density of free massless 
Wilson fermions vs. N4 for three aspect ratios in 
figure ^ As in the discretized case of section ^, 
approaching the A'4 —^ 00 limit at fixed small 
aspect ratio underestimates the SB result. How- 
ever, the most obvious feature is the large over- 
estimation of the SB result at smaller A'4. At 
the smallest A'4, at least, we now can understand 
that the energy density is cutoff by the scarcity 
of lattice momenta. 

We propose a simple model to reveal the source 
of the overestimate for Wilson fermions. Without 
the Wilson term, the energy density of free naive 
lattice fermions would be 16 times larger than the 
SB result. The Wilson term gives 15 doublers a 
mass proportional to the cutoff. But from fig- 
ure |l|, we expect these now massless doublers to 
make a significant contribution to the overall Wil- 
son energy density until A'4 > 8. To emphasize 
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Figure 4. Discretized energy density with fifteen 
"doubters." 

the point, in figure ^ we compare the free Wilson 
energy density with that of one massless and fif- 
teen massive flavors using the discretized energy 
density (|). While more effort could be made 
to adjust the masses of the discretized doubters 
to make the two curves agree, from the common 
shape we believe this is a good explanation of the 
Wilson excess. 

4. Anisotropic domain wait action 



The energy density is formally given by 
1 dlnZ 



— £ = 



V c»r-i 



(3) 



V 



In order to vary the temperature continuously 
while holding the volume fixed to compute the 
derivative, we must allow the lattice spacing in 
the temporal direction, a4, to differ from the spa- 
tial direction, a. The ratio of the spacings is the 
anisotropy = a/a^. Of course, once the deriva- 
tive has been computed we can take the isotropic 
limit, ^ — > 1. 

To derive the form of the anisotropic domain 
wall and overlap fermion matrix, we choose the 
interpretation of the extra dimension in the do- 
main wall formalism as a flavor space. We start 
from Klassen's anisotropic Wilson action with- 
out the anisotropic clover terms as the domain 



wall action is 0{a) improved in the Ls — > 00 limit. 
Using the conventions of the Columbia group jj] 
the anisotropic domain wall action is 



{x, x') = E y - U,{x)5.,+ -,.,., 

+ [mo-r(3eV + ei'4)] (5) 

where the are the bare velocities of light with 
vi = V2 — = the anisotropy is f = aj a/^^ and 
the exponents are set Xo z\ — zi = zj, = z = ^ 
and Z4 — 1 although other choices are possible if 
they satisfy the constraint Z4 — z = 1 [|| . The fla- 
vor mixing matrix D-^{s, s') is the same as written 
by the Columbia group [Q without any extra co- 
efficients since the normalization of the fermion 
fields was shifted so that no coefficients appear 
with the mo term in the action. The domain 
wall action includes a Pauli-Villars subtraction at 
TTif — 1 to cancel the bulk infinity from the Ls 
heavy flavors at the mass scale rriQ. 

Following Neuberger the anisotropic form 
of the overlap action can be derived as the Ls — + 
00 limit of the anisotropic domain wall action. 
However, the form of the action is easy to guess: 
just replace the isotropic Wilson matrix with the 
anisotropic one in the overlap Hamiltonian. 

5. Domain wall/overlap energy density 

From (H), the energy density is computed ex- 
plicitly using 



ji4 



- Tr 



D 



(6) 



Our approach was to compute the derivative 
£,d/d£,Dp analytically and then diagonalize (|^) 
into blocks indexed by four-momentum p^^ us- 
ing Fourier transform. Each block matrix Dp{p) 
was inverted numerically and then the traces were 
summed over all four-momenta. In general, (||) is 
computed twice, once on a (CA^4)^ X Na lattice 
and again on a {(^N^Y lattice for the T = vac- 
uum subtraction. Anti-periodic boundary condi- 
tions were imposed in the temporal direction in 
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Figure 5. Domain wall and overlap energy den- 
sity at mo = 1. 

both cases. For domain wall fermions, the cal- 
culation is repeated with the Pauli-Villars mass 
and the difference is used, schematically edwf — 
[epiT) - £f (0)] - [spviT) - epviO)]. 

The energy density for domain wall fermions at 
various Lg and overlap fermions are shown in fig- 
ure ||. The choice mo = 1 is the optimal free 
field value in the sense that the surface mode 
with = is completely localized to the do- 
main wall. However, surface modes with other 
momenta will have finite exponential decay rates 
so the Ls dependence is due to these modes. The 
energy density at Ls = 10 is nearly converged on 
the Ls oo limit. Any small differences between 
the domain wall and overlap at smaller A'4 may 
be due to the different treatment of the spacing 
of sites in the extra dimension. This is currently 
under study. 

The Wilson energy density is also shown in fig- 
ure ||. Based on arguments of section ||, the re- 
duced excess energy density for domain wall and 
overlap fermions indicates the participation of the 
heavy doublers in the thermodynamics is sup- 
pressed due to Pauli-Villars term. In figure ^ we 
show the energy density of domain wall fermions 
before Pauli-Villars subtraction and separately 
the energy density of the Pauli-Villars regulators. 
Before subtraction, the excess energy density due 
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Figure 6. Domain wall fcrmion and Pauli-Villars 
energy densities. 

to heavy modes is essentially Ls times larger than 
the Wilson case, again shown as a dashed line. 
That the domain wall density of order one appear 
as the difference of two densities of order Ls sug- 
gests tinkering with the Pauli-Villars subtraction 
may reduce heavy contributions further. How- 
ever, it is unclear how worthwhile a lengthy study 
of possible Pauli-Villars terms would be when the 
benefits of such a finely tuned action may vanish 
in the interacting case. 

6. "Doubly" regularized fermions 

As pointed out by Neuberger lattice 
fermion actions with Pauli-Villars subtractions 
are "doubly" regularized, once by the lattice 
spacing and again by the Pauli-Villars pseud- 
ofermions. In section ||, we claim the Pauli-Villars 
suppression of heavy contributions to the domain 
wall energy density works so well the end result 
is an improvement over Wilson fermions quite 
apart from any issues of chiral symmetry. How- 
ever, there is no reason we can't attempt to use a 
Pauli-Villars regulator to achieve the same result 
for Wilson fermions. 

In figure ^ we show the Wilson energy density 
before Pauli-Villars subtraction, as has already 
appeared in the last several figures as open cir- 
cles, the energy density of the proposed Pauli- 



^=4 (aspect ratio) 



■— ew(0)-ew(1) 




N, = (aT)- 



Figure 7. Energy density of "doubly" regularized 
Wilson fermions. 

Villars regulators with a mass of mo = 1 as open 
diamonds and the doubly regularized Wilson en- 
ergy density appears as filled squares. After the 
subtraction, the contribution of heavy modes to 
energy density is substantially reduced. Fine tun- 
ing of the Pauli- Villars mass to make the energy 
density agree with the Stefan-Boltzmann limit at 
a given N4 is certainly possible, but it is again 
unclear how to do this in the interacting case. 

Another potential advantage of doubly regu- 
larized fermion formulations is that the quenched 
theory is recovered when the fermion mass is set 
to the Pauli- Villars mass. An example of its use- 
fulness is a proposed calculation of the dynamical 
QCD equation of state in the deconfined phase us- 
ing the integral method. As the integration con- 
tour originates in the strong coupling regime and 
crosses through the transition region, an inordi- 
nate amount of work must be done before any- 
thing is known of the deconfined region. As the 
quenched theory appears as a different contour 
of fixed finite mass, one could use the integration 
contour shown in figure ^. As the quenched equa- 
tion of state is already well known for several A'4 
(2I , all dynamical simulations are performed only 
at weaker couplings where fermion algorithms are 
usually better conditioned. 



quenched {m,=mpy) 




/ 

/ 


7 




/ 

/ ' 


/ 

/ 


1 

1 

1 

1 




/ 

/ 

/ 






/ 


' /' 

/ 




/ 

/ 

/ 4 

/ 


/ 

/ 

/ 









Figure 8. Integration contour for deconfined 
equation of state 

One subtlety of integrating "off the quenched 
contour" is the importance of computing the T = 
vacuum subtraction at as low a temperature as 
possible. As before, finite temperature observ- 
ables are computed on x A^4 lattices and 

r = subtractions are computed on (C^4)'* lat- 
tices. On the quenched contour (aT)^^ — ^iV4 
lattices may be confined, but while integrating 
down at fixed coupling, one may enter the (^A^4 
transition region. To avoid this, a stair step con- 
tour could be used. 

We would like to thank N. Christ, R. Mawhin- 
ney, T. Klassen and P. Vranas for useful discus- 
sions in the early stages of this work. 
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